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ABSTRACT 

We consider the passage of a beam of non-thermal electrons through 
the flaring solar atmosphere, paying particular attention to the require- 
ment that the beam be stable to the generation of plasma turbulence. We 
then compute the ratio of energy losses due to reverse current ohmic 
heating, and heating by Coulomb collisions, respectively, for the greatest 
flux which can pass stably through the atmosphere. We show that this 
ratio is determined by the low energy cutoff of the beam, by the electron 
temperature of the ambient atmosphere, and by the electron to ion tempera- 
ture ratio 9 . It is also independent of the atmospheric density. The 
results show that ohmic energy losses are undoubtedly important in the 
initial transient state, in agreement with other authors, but that their 
role is debatable in the flare atmosphere, depending on the value of 0 
appropriate. Expected values for 9 during the impulsive phase of the 
flare indicate that reverse current ohmic energy losses, and their con- 
sequent effects on the electron beam dynamics and the hard X-ray 
bremsstrahlung emission, may not be as important as previously suggested; 
however, a fully time-dependent analysis of the beam-target interaction 
is necessary to fully resolve the issue. 

Subject headings : hydromagnetics - particle acceleration - Sun: corona 


Sun: flares - Sun: X-rays 
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I . INTRODUCTION 

The necessity for a beam-neutralizing reverse current in a thick 
target electron-heated scenario of solar flares is by now well established 
(e.g. Hoyng, Brown, and van Beek 1976; Brown and Melrose 1977; Knight 
and Sturrock 1977; Hoyng, Knight, and Spicer 1978; Emslie 1980). In addi- 
tion, detailed knowledge of the various processes affecting the passage of 
an electron beam through the solar atmosphere, and their relative importance, 
is essential in order to infer the characteristics of the accelerated elec- 
tron population (and so place constraints on the flare primary energy 
release process) from radiation signatures such as hard X-ray emission. 

For this reason it is important to ascertain as accurately as possible the 
effect of the reverse current on the evolution of the distribution of acce- 
lerated non-thermal electrons. 

A quantitative treatment of the ohmic heating of the flare atmosphere 
resulting from the passage of the reverse current through the resistive 
ambient plasma was performed by Emslie (1980). He showed that these ohmic 
losses can, under certain conditions, result in considerable modification 
to the dynamics of the electron beam (compared to those calculated under a 
purely collisional treatment - see Emslie 1978). In addition, he showed how 
the features (e.g., intensity, spectrum) of the hard X-ray bremsst rah lung 
produced by collisions of the beam electrons on ambient protons may also 
be significantly affected. In his analysis, however, he assumed, without 
rigorous justification, that the reverse current drift velocity was such 
that plasma instabilities did not develop. Hoyng, Knight, and Spicer (1978) 
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did consider the question of the stability of the reverse current (against 
the generation of ion-acoustic turbulence) and concluded that if a flux 
F unstab non_t ^ ern,a ^ electrons is injected with an associated reverse 
current which would be unstable, then the instability acts to reduce this 
injected flux to a marginally stable value F cr i t < F uns tab ^ see a * 80 
Brown and Melrose 1977; Manheimer 1977; Spicer 1977). 

It is therefore of interest to consider the relative roles of ohmic 
(i.e., due to the passage of the [lov ocity] reverse current through the 
finite resistivity background plasma) and collisional (i.e., due to the 
interactions of the [high velocity] beam electrons with the ambient particles) 
heating for a beam of non-thermal electrons whose flux, at injection, is 
such that the beam is marginally stable to the generation of plasma turbu- 
lence. We find (ill) that the ratio of the energy loss rates by both pro- 
cesses is dependent on the low energy cutoff in .the injected electron spec- 
trum, on the electron temperature of the target, and on 3 , the ratio of 

electron (T g ) to ion (T^) temperatures in the target. The first two of 

these parameters are quite well established (to within factors of two); 
however, the value of 6 is not so well established, and the results 
are found to be quite sensitive to this parameter. We thus find (S III) 
that, ohmic energy losses will only be important if 0 is comparatively 
low (< 3), while for larger values of 0 a purely collisional treatment 

'Sj 

is adequate. This threshold value of 0 is then compared with those 
expected under a variety of conditions, allowing us to draw qualitative 
conclusions regarding the role of reverse current ohmic heating at 
various stages in the flare. 

or ' ' ' • - 

OF FoO.H Q./.LiTY 
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II. THE RATIO OF OHMIC TO COLL1SIONAL ENERGY LOSSES 
FOR A MARGINALLY STABLE ELECTRON FLUX 

We shall consider Che injection of a beam of non-thermal electrons 
with energy spectrum 

F / E 

W * ( 6_1 ) e^(e^) electrons cm“ 2 s“ 1 keV~ 1 , ^ 

for E q > E^ , the low energy cutoff (keV) for the beam. (Note that 
is thus the total injected flux.) Ohmic energy losses are greatest near 
the acceleration site (see Figure 3 of Emslie 1980) and so we will com- 
pare ohmic and collisional energy loss rates at that point. Using 
equation (34) of Emslie (1978) to calculate the collisional losses 
(assuming a fully ionized target) and setting the reverse current losses 
equal to nj 2 ■ ne 2 F^ (n * plasma resistivity, j * reverse current 
density), we find that the ratio of ohmic to collisional losses is given 
by 

6 - — , ( 2 ) 

2ir e 2 (6-l)An 

-3 

where e is the electronic charge (e.s.u.), n the target density (cm ), 
and A the Coulomb logarithm (Spitzer 1962). 

Since the electron flux produced by any source will be limited to the 
marginally stable value, as discussed by Hoyng, Knight, and Spicer (1978; 
see discussion in §1), we may obtain an upper limit to B , B q . by 
setting F^ ■ F 9 at mar 8 inal stability. Further, since 

outside the electron source, i.e. in the region of beam propagation, the 
flux will always be lower than this (due to collisional and reverse current 
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ohmic attenuation processes [see Emslie 1980]), we may set all plasma trans 
port coefficients to their classical values (Spitzer 1962). In particular, 
therefore, the resistivity n will assume the form 


n 


3/2 1/22 
it m e A 

e 


2(2kT } 


3/2 


(3) 


where m is the electron mass and k is Boltzmann’s constant, F . 

e crit 

is given, for marginal stability, by 


F 


crit 


anv , 
e 9 


(4) 


1/2 

where v ■ (kT /m ) is the electron thermal velocity and a is a 
e e e 

function of the electron to ion temperature ratio 0 • Substitution of 
equations (3) and (4) into equation (2) yields an upper limit to the 
contribution of ohmic losses (relative to collisional ones) in the target: 


6 < 6 0 



(5) 


Note that this expression is independent of the background density n . 

The form of a(9) depends on the conditions in the region where the 
electrons stream. Duijveman, Hoyng, and Ionson (1981) have 

considered the effect of both ion-cyclotron and ion-acoustic wave genera- 
tion on a driven electron flux. They find that for 0 > 8 the marginally 
stable flux is determined by the ion-acoustic turbulence threshold; the 
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relevant form is then as in Figure 8 of Fried and Gould (1961) 1 . For 

1 Some authors (e.g., Brown, Hayward, and Spicer 1981) have used an analytic 

3/2 

approximation for a(0) , which scales as (1 + A0 exp [-0/2]), where 
A • (Op/mJ 1 ^ 2 e~ 3 ^ 2 , m p being the proton mass. This function in fact 
fits the exact numerical results of Fried and Gould (1961) only for 0 > 6 . 
It also exhibits a maximum turning value of * 0.28 at 0 ■ 3; there is 
clearly no justification for use of this maximum value, especially when one 
considers that it lies w^ll outside the range of applicability of the 
approximation. 


0 < 8 , however, the situation is somewhat more involved. For such low 
values of 0 the onset of ion-cyclotron turbulence precedes the onset 
of ion-acoustic turbulence (i.e. the critical drift velocity for the 
former is smaller; see Figure 1A of Duijveman, Hoyng, and Ionson 1981); 
which stability threshold is relevant then depends on the turbulence level 
at which the ion-cyclotron waves saturate. If the ion-cyclotron wave 
turbulence does saturate (i.e. the saturation level is sufficiently low), 
then one may drive electrons to velocities larger than the ion-cyclotron 
stability threshold; a is then determined by the ion-acoustic stability 
threshold, as for 0 > 8 . If, on the other hand, the waves do not 
saturate (due to a high saturation level), tnen the ion-cyclotron stability 
threshold determines a . Since the saturation level of ion-cyclotron 
waves in solar conditions is not well known (Duijveman, Hoyng, and Ionson 
1981), it is in practice uncertain which of these two possibilities for 
0 < 8 prevails. 
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Figure 1 


In Figure 1 we show the dependence of 6 Q on 0 , determined from 
equation (5) and the above considerations of the value of o appropriate, 
for a typical 6 ■ 4 and for various values of Ej/T.^ , where E^ is 

measured in keV and T7 is the coronal electron temperature in units of 

10 7 K . The solid lines correspond to the ion-acoustic threshold and 

the dashed lines to the (possibly irrelevant, depending on the turbulence 

saturation level) ion-cyclotron threshold. Note that when 0 < 8 and the 

ion-cyclotron wave saturation level is low, the resistivity 'l^, associated 

with the saturated ion-cyclotron turbulence, is less than the Spitzer (1962) 

classical value (see Figure 5 of Duijveman, Hoyng, and Ionson 1981). Thus 

classical resistivity remains a valid approximation, so that equation (5), and 

also the solid curve in Figure 1, still give a good estimate of 6 under 

o 


these conditions. 
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III. DISCUSSION 

In order Co interpret Che resulcs of Figure 1, we muse firsc assign 
an appropriace value eo Ej/T 7 . T 7 Increases from * 0.3 in che preflare 
stace (Noyes 1971) Co a maximum of around (2 ± 1) (e.g. Dadowe, Hudson, 
and Pecerson 1974). Power-law X-ray ijpectra have been observed down Co 
c • S keV (Kane eC al. 1979); when one cakes inco accounc che chermal 
concribucion co Che phoCon flux ac chese energies, one finds ChaC chis 
corresponds co Ej * 7 - 8 keV (lrown, Hayward, and Spicer 1981). He 
shall Cherefore adopt Ej/T 7 * 25 for Che preflare scace, and Ej/T 7 * 10 
as an upper limic for Che flare lcself. 

Considering firsc che preflare acmosphere, reference co Figure 1 
shows ChaC 6 Q exceeds unicy for all values of 8 < 5 , whacever che 

marginally scable scace in che cargec may be. Thus, 

since 6 s 1 in Che preflare scace, reverse currenc ohmic energy losses can 

easily dominace over collisional losses, confirming che conclusions of Knighc 

and Sturrock (1977). However, in the flare atmosphere, S Q exceeds unity 

only for 0 < 3 in the case of an ion-acoustic turbulent threshold, and 

2 

only for 8*1 for an ion-cyclotrou turbulent threshold. Because 

2 The resulcs of Emslie (1980) , showing ChaC reverse current effects are 
strong for large injected fluxes, tacitly assume that such large fluxes 
are scable to the generation of plasma turbulence. T lis implies that a 
reverse current with a large drift velocity be permitted to pass stably, 
which in turn implies that 8 be close to unity (Fried and Gould 1961; 
Duijveman, Hoyng, and Ionson 1981). 
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Ccuiosb co'.lisicr.* and ohmic heating both preferentially heat electrons 
opposed to ions, it would appear that, since the initial ion tempera- 
ture is only of order 10 6 K, the appropriate value of 0 is much larger 
than either of these two values, implying 6 0 << 1 . 

However, the heading of the plasma by both collisional and ohmic 
processes is counteracted by conductive and radiative cooling of the 
heated region. If a steady state situation obtains, then the external 
source terms in the equations controlling the behavior of T £ and 
will vanish and T # and will be governed simply by the equili- 

bration equations 



where 

- 12.6 T 3 ' 2 tf 1 seconds 

e 


( 6 ) 

(7) 


is the (classical; see §11) electron-ion temperature equilibration time 

(Spitzer 1962). The solution of equations (6) and (7) has been obtained by 

Sivukhin (1966, his equation [9.17] and Figure 6); for an Initial state in 

which 6 is large (due to preferential heating of electrons; sec remarks 

above) one finds that 6 falls to a value £ 3 in a time r* * 0.3 T and 

o 

to a vftluft a 1 in T* * Tq , where is the initiei temperature equilibra- 

tion timescale. The first of these t* values corresponds to Bq ■ 1 

when the ion-acoustic turbulent threshold applies and the latter to 6 ■ 1 

o 

when the ion-cyclotron turbulent threshold applies (Figure 1). 

We thus see that if the beam injection time is long comparer 1 to 
t* and to th* time taken to achieve a balance between source and sink 
terms in the energy equation, then a state of equal electron and ion 
temperatures will be reached, and reverse current ohmic energy losses 


' ‘ - r 

Q 


i 

i V 
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can in fact become important again, at least near the electron injection 
site (Figure 1; Figure 3 of Emslle 1980). If, on the other hand, either 
t * ur the tine to achieve local energy balance is much larger than the 
duration of the electron inp«.t, then 0 will remain large and so reverse 
current ohmic energy losses will be unimportant (relative to collisional 
ones) everywhere in the target. 

These considerations reveal that the comparative roles of reverse 

current ohmic, and direct collisional, energy losses are determined by 

the target density, not (as might be expected) directly through its 

appearance in the collisional energy deposition rate (erg cm ^ s~*) , but 

through its effect on the local energy balance in the flaring corona and 

on the electron-ion temperature equilibration time t (eq. [7]). Typical 

preflare coronal densities are > ..O 10 ca- so that (assuming hydro- 

dynamic effects are small) t < 40 a. Ve thus see that t* could be 

o 

as low as 10s (when the icn-acoustic turbulent threshold applies); this 
is comparable to the duration of an "Elementary Flare Burst" (de Jager 
and de Jonge 1978), allowing significant temperature equilibration to 
occur should balance between heating and cooling ter-is be achieved within 
such a timescale. If, however, a state of local energy balance takes 

longer than this to be attained, then 8 will remain >> 1 and so 8„ 

0 

will remain << 1 , as discussed above. 

In summary, therefore, the relative roles of reverse current ohmic, 
and collisional heating are determined by the sagnitude of the injected 
non-thermal electron flux, in addition to the parameters of the target 
atmosphere. An upper limit to this injected flux is that which is 
marginally stable to the generation of plasma turbulence, either ion- 
acoustic or ion-cyclotron as appropriate. In a preflare atmosphere, an 


11 


injected flux of thie magnitude produces strong ohmic heating, in - s ree- 
ment with the conclusions of Knight and Sturrock (1977). T:. a flare 

atmosphere, however, the importance of ohmic energy locses is determined 
principally by the critical reverse current drift velocity and hence by 
the electron to ion temperature ratio 6 . This in turn depends on the 
relative sizes of three timescales: the elr tron beam lifetime, and 

the timescales for achieving local energy balance, and equilibration of 
electron and ion temperatures, respectively. If a state of local energy 
balance is reached while the electrons are still being injected, then 6 
will approach unity within the order of an electron-ion temperature 
equilibration time, after which reverse current ohmic energy losses will 
be an important consideration, at least near the electron injection point, 
for suitably large (but stable) injected electron fluxes (Emslie 1980). 

On the other hand, if a state of local energy balance cannot be reached 
this quickly, then 6 will remain large, so rendering reverse current 
ohmic energy losses negligible compared to colllsional losses throughout 
the entire atmosphere for any stable electron flux. A detailed study 
of the time-dependent response of the solar corona to a non- thermal 
electron energy input is necessary in order to resolve these questions. 

In addition, it is important «.o assess whether different "Elementary Flare 
Bursts" (corresponding to different "spikes" in the hard X-ray flux-versus- 
time profile -- see da Jager and de Jonge 1978) result from repeated re- 
energization of the same region, or to successive energizations of differ- 
ent regions (see Karpen, Crannell, and Frost 1979; Emslie 1981). In the 
former case, the time available for electron-ion temperature equilibration 
will be relatively large, while in the latter case each burst must be 
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considered individually, with preflare initial conditions, thus requiring 
that the electron-ion temperature equilibration time be less than the 
duration of the Elementary Flare Burst if reverse current ohmic losses are 
to become an important energetic consideration at some stage in the burst. 
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FIGURE CAPTION 


Fig, 1 - Ratio of reverse current ohmic heating to collisional energy 
deposition, 8 q , as a function of the electron-ion temperature ratio 
9 . The values of 8 q have been calculated for a non-thermal electron 
flux whose associated reverse current, at the point of injection, is 
marginally stable to the generation of plasma turbulence. They are 
also evaluated at this injection point, and so represent an upper limit 
to the importance of reverse current ohmic heating in the atmospheric 
energy balance (Emslie 1980). The dashed curves correspond to the 
marginally stable threshold for ion-cyclotron turbulence, and the solid 
curves to the marginally stable threshold for ion-acoustic turbulence, 
applicable either when 9 > 8 or when 0 < 8 and the ion-cyclotron 
wave saturation level is low (Duijveman, Hoyng, and Ionson 1981). In 
both cases 8 q depends linearly on the ratio of the low energy cutoff 
in the electron spectrum (keV) and the electron temperature of the 
target atmosphere (T^ in units of 10 ^K), as shown. 
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